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ABSTRACT 

The recently reported microlensing events of the LMC have caused much excitement, and 
have been interpreted as due to 'dark objects' (MACHOs) in the halo of our Galaxy. It 
is shown here that the stars within the LMC play a dominant role as gravitational lenses 
and can indeed account for the observed events. 

For observations within the bar of the LMC, the probability of microlensing being caused 
by a star within the LMC is found to be ~ 5 x 10~^. Outside of the bar, the probability 
of microlensing being caused by a star in the LMC is 4 to 12 times lower. The MACHO 
event (Alcock et al. 1993) and one of the EROS events (Aubourg et al. 1993) lie within 
the bar for which the probability of microlensing is consistent with being caused by an 
object within the LMC. 

If the stars within the LMC play a dominant role as lenses, the events should be concen- 
trated towards the center of the LMC. On the other hand, if MACHOs play a dominant 
role as lenses then, for a given number of monitored stars, the events should be uniformly 
distributed over the whole area of the LMC. Thus the galactic and the LMC lenses can be 
statistically distinguished in most cases. 

It is further shown that, under certain conditions, the light curve of an event caused by a 
star within the LMC would be different from the one caused by a MACHO. This can also 
be a distinguishing signature, and seems to have been observed in case of the MACHO 
event. The fit discrepancy near the peak which the authors say "is not yet understood" 
would be a natural consequence if the event is caused by an object within the LMC, which 
further indicates that the lensing is due to a low mass star within the LMC itself. 
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1. Introduction 



Most galaxies are known to have flat rotation curves (Begeman, 1987; Sanders and Bege- 
man 1994). In the outer parts of the galaxies in particular, the visible matter falls far short 
of what is required to explain the flat rotation curves. In our own Galaxy, the rotation 
curve is observed to be flat up to at least 16 kpc from the center (Fich et al. 1989). The 
scenarios proposed to explain the flat rotation curves are either departure from Newtonian 
dynamics in large scale, or presence of dark matter (for a review see Sanders 1990 and 
Broeils 1991). In the latter and more conventional scenario, the rotation curves would 
imply that a significant part of the mass of the galaxies, including our own, resides in the 
halo in some form of dark matter. The nature of this dark matter has been hypothesized 
to be in either of the two forms: MACHOs (i.e. massive compact halo objects which is 
a collective term for 'Jupiters', brown dwarfs, red dwarfs, white dwarfs, neutron stars or 
black holes); or elementary particles such as massive neutrinos, axions or WIMPs (weakly 
interacting massive particles such as photinos etc.). One important distinguishing fac- 
tor of the MACHOs is that they can have detectable gravitational effects and can cause 
"gravitational microlensing" of background stars. 

Einstein (1936) was the first to point out that a star can act as a gravitational lens for 
another background star, if the two are sufficiently close to each other in the line of sight. 
Given the observational capabilities of that time, Einstein had however considered this to 
be a purely theoretical exercise and had remarked that there was "no hope of observing 
such a phenomenon directly" . 

The observational capabilities have improved remarkably since then. More recently, 
Paczyhski (1986) worked out the probability of such microlensing events by MACHOs 
and showed that if the halo of our Galaxy is made up of MACHOs, the probability of 
finding them through microlensing is 5 x 10~^, independent of their mass distribution. He 
suggested an experiment to look for such events using the LMC stars. Further details of 
this idea were worked out by Griest (1991). Such an experiment was taken up by 2 groups 
who have reported their first results (Alcock et al. 1993; Aubourg et al. 1993). 

2. The experimental results 

The results from the MACHO experiment (Alcock et al. 1993) are from the monitoring of 
1.8 million stars in the bar of the LMC, both in blue and red wavelengths, for a period of 
one year, which led to the detection of one event. More recently, it has been reported that 
the analysis of all the data, which involves monitoring of 8.3 million stars in the region 
of the bar over 2 years, has led to the detection of 4 events (Cook et al. 1994). The 
EROS experiment (Aubourg et al. 1993) involves 2 programs: one involves monitoring of 
8 million stars over a field of 5 x 5 degrees using Schmidt plates in red and blue, with a 
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sampling rate of one or two per night over a period of 3 years. The second program involves 
monitoring of 100.000 stars in red and blue wavelengths, over an area of 1 x 0.4 square 
degrees with a sampling rate of about half an hour, for about 7 months. The analysis of 
40% of data of the first program led to the detection of 2 events, and no event was detected 
in the second program. The details of the published events are listed in Table 1. 

3. Concerns if the microlensing is caused by MACHOs 

If these are genuine microlensing events by galactic halo objects, this could potentially 
herald an end to the problem of dark matter. But this has also been the cause for some 
concern: if these events are caused by the objects in the halo, it would lead to some 
problems in terms of stellar evolution and the theory of galaxy formation (Hogan, 1993). 

Recent work on deuterium abundance derived from the observations of a quasar suggests 
that the dark halos of galaxies may well be nonbaryonic (Songaila et al. 1994). 

In case of the MACHO event, which has the highest signal to noise ratio, there is a 
discrepancy between the best fit and the observed light curve near the peak of the light 
curve. This consistently repeats both in the blue and in the red, and according to the 
authors this "fit discrepancy near the peak is not yet understood" . Furthermore, the rate 
of microlensing events as observed seems to be lower than expected from the dark halo, and 
it is important to estimate the rate expected from the stars that are known to exist in our 
galaxy and in the LMC. The analyses for the galactic stars and the LMC halo have been 
done, but the importance of the LMC stars had been overlooked so far. Sahu (1994) has 
recently shown that the stars within the LMC indeed play a dominant role as gravitational 
lenses, the full details of the probability analysis is presented here. It should be noted here 
that, for the microlensing events of the galactic bulge stars (Udalski et al. 1993; Udalski et 
al. 1994), the stars within the bulge also play an important role as gravitational lenses, the 
analysis of which has recently been done by Kiraga and Paczyhski (1994). Furthermore, 
this paper presents the analysis of the light curve and shows that if a star within the LMC 
acts as a lens, the source can be extended which reproduces the light curve better. Before 
attempting to analyse the light curve in detail, let us first see in detail the probability of 
this scenario. 

4. Analysis of the probability 
4.1 Mass of the bar and stellar density 

To calculate the probability of the microlensing being caused by a star in the LMC itself, 
we need to know the stellar mass density in the LMC. Let us first confine ourselves to the 
bar of the LMC. From the surface luminosity, it is estimated that the observed luminosity 
of the bar is about 10 to 12 % of the total observed luminosity of the LMC in the optical 
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wavelengths (De Vaucouleurs and Freeman 1973, Bothum and Thompson 1988). To see 
whether the effect of extinction is important and to account for it, we note that the 
extinction in the region of the bar, as derived from the 100 /im. IRAS maps is between 12 
and 24 MJy per steradian with an average value of about 17 MJy per steradian (Schwering 
and Israel 1990), which translates to an average extinction of about 1.5 mag in the V 
band (Laureijs, 1989). This value is consistent with the value deduced by Hodge (1991) 
from background galaxy counts and is also consistent with the value derived from the 
observations of the most reddened stars in the region (Isserstedt 1975). The average 
extinction in the outer regions is about 0.3 to 0.4 magnitudes. It is perhaps worth noting 
here that, as we will see later, errors in extinction measurements have only very little effect 
on the final probability. To see how the extinction affects the mass, let us assume that the 
extinction is uniform in depth, and let d be the total depth. Thus for any line of sight, if 
A^; is the total extinction, the ratio of the observed to the true luminosity at that point 



1 1 ,,-n.oi6A^ 



Ltrue d Jq 0.916 a 

where z is the direction of the line of sight. 

Using Ay=1.5 mag for the bar, and 0.3 to 0.44 mag for the region outside, we get 
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Now, using the fact that {Lobs)i,ar/ iLobs)tot ~ ^-^ 0.12, and after a little algebra we get, 

{Ltrue)bar ^ q -^^^ q -^^ 
[LtrueJTot 

Thus the true luminosity of the bar is about 14 to 18 % of the total luminosity of the 
LMC. The total mass of the LMC, as determined from various methods, ranges from 6 x 
10^ to 15 X 10^ M© (De Vaucouleurs and Freeman, 1973). (I n case of the LMC, there is 
no significant discrepancy between the mass derived from the observed luminosity and the 
mass derived from the velocity dispersion.) Assuming that the mass to light ratio in the 
bar and the outer parts are the same, we get the mass of the bar to be in the range 1 x 
10^ M© to 2.7 X 10^ M©. In our subsequent calculations, we will round off the mean and 
use 2 X 10^ M© as the mass of the bar. Considering the fact that the LMC is gas poor 
and only about 5% of the LMC mass is neutral hydrogen ( Wester lund 1990, Israel and de 
Graauw 1991, Rohlfs et al. 1984), we will neglect the contribution of gas to the mass and 
assume this entire mass to be made up of stars. 

4.2 Expected probability of microlensing in the region of the bar 
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To calculate the probability of microlensing, let us assume that there are N^o^ stars being 
monitored. We note that, the stars that are farther away, have more stars in front of them, 
thus having a larger probability of being microlensed. But the observed number density 
of stars becomes less at larger depths because of extinction. Thus, in order to correctly 
calculate the probability, we must know how the extinction affects the observed number of 
stars with depth. To make life easier, we can take the current monitoring programs and 
the magnitude limits to calculate the probability. We note that the limiting magnitude 
of the current surveys are about 20 to 21, which at the distance of LMC, corresponds to 
an absolute magnitude of 1.5 to 2.5. So the magnitude of the stars that can be observed 
at the near side is 1.5 to 2.5 whereas the magnitude of the stars at the far end is 1.5 
magnitude brighter. If the distribution of stars among different spectral types is assumed 
to be similar to what is observed in the solar neighborhood, then the difference in stellar 
number density from front end to the back is about 3 (Allen 1973). Similar results are 
found from HST observations of LMC clusters (Gilmozzi et al. 1994). (If the extinction 
is 3 magnitudes, then this value is about 10. As it turns out, the effect of extinction is 
small. If the extinction is larger, the mass of the bar becomes larger thus increasing the 
probability of microlensing. But the extinction effectively reduces the observed number of 
stars which are at a larger depth which have the larger probability of being microlensed, 
thus decreasing the probability. These 2 effects compensate for each other the net effect 
is that the probability decreases with extinction, but only slightly. If the extinction is 3 
magnitudes instead of 1.5, the net probability changes only by less than 50%. Furthermore, 
at the magnitudes considered here, the age is not expected to play a significant role, so to 
a good approximation, the spectral distribution of the stars can be assumed to be similar 
to that in the solar neighborhood for the probability calculation.) 

Assuming the extinction to be uniform in depth, we can express the observed number of 
stars at any layer dz, at a depth of z, as 



(2) 



where A^^ is the observed number of stars per unit depth in absence of extinction i.e, the 
observed number of stars per unit depth at the near side. If Ntot is the total number of 
stars being monitored, then 



Ntot = / N,3-^^dz = N,d 
Jo 

Substituting eq. 3 in eq. 2, we get 
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The fraction of area covered by the Einstein rings af all the individual stars lying in the 
front of this layer can be expressed as 

Af{z) = [ TTR%{z)n{z)dz (5) 
Jo 

where n{z) is the stellar number density at depth z. Since n = p/m and R% oc m where p 
is the stellar mass density and m is the mass of the star 

c Jo 

which is thus independent of the mass distribution. Assuming the stellar density to be 
uniform with depth 

27TGpZ^ 

Fom Eq. 4 and 7, the instantaneous probability of observing one microlensing event, when 
Ntot stars are being observed, is 

P = jyo,s{z)Af{z)dz = ^yi^/'s-t^^d^ (8) 

Strictly speaking, this is valid for P « 1, which is the case here. The optical depth 
r, which is equivalent to the probability of any star being microlensed at a given time is 
P/Nioti is easily obtained by solving eq. 8 as 

/ (9) 



For the value of p, we can substitute, 

M 



(10) 



^ LWd 

where L, W, d and M are the length, width, depth and mass of the bar respectively. 

Assuming L=3000pc, W=d=600pc and M = 2 x 10^ M0, and substituting eq.lO in eq.9, 
we get, 

nar ^ 5 X 10"^ (11) 

with some uncertainty due to the uncertainties in the size and mass of the bar. It must be 
mentioned here that the value of the optical depth calculated above is an average value for 
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the whole bar. Since the luminosity gradient in the bar itself is about 30% from the center 
to the outer parts (Bothum and Thompson 1988), the optical depth would be slightly- 
higher in the central region of the bar, and slightly lower in the outer region. 

It is clear from Eqs. 9 and 10 that t cc d. In absence of any more accurate information, 
we have used the depth to be the same as the width and have neglected any inclination 
effect. In the literature one finds that the bar can be as thin as the disk itself (Binney 
and Tremaine 1987). So, considering the fact that the thickness of the LMC disk is about 
300pc, and the inclination is about 45 degrees (Westerlund 1991), the depth can be about 
450pc in the extreme case, thus the effect due to the uncertainty in d can be 25% at the 
most. 

Carrying out an identical analysis for the region outside the bar it is easy to see that, if 
the depth of the LMC disk is between 100 to SOOpc (Feast, 1989), the optical depth in the 
region outside the bar is about 4 to 12 times smaller. 

4.3. Comparison with observed probabilities 

The MACHO event lies in the central region of the bar. Keeping in mind the uncertainties 
involved in the estimation of optical depth on the basis of the single published MACHO 
event, we can only say that it seems to be well below the optical depth of 5 x 10~^ 
expected from a dark halo made up entirely of MACHOs, and is consistent with the 
optical depth calculated above. Extending an identical analysis for the more recent results 
of 4 events reported from the monitoring of 8.3 million stars over a period of 2 years, the 
estimated observed probability is found to be indeed very close to the value calculated 
above, although there is some uncertainty in the value of the detection efficiency. In case 
of the EROS events, one event lies in the outer region of the bar, the other is far from 
the bar, and the optical depth has been estimated to be higher (~ 2 x 10~^, Gould et al. 
1994). But Gould et al. (1994) also suggest that there may be some systematic effects and 
the detection efficiencies may be uncertain. The situation will be clearer as more events 
are observed, particularly events with higher magnifications (see below), and we must wait 
for more events to be observed before making any direct comparison with the calculated 
optical depth. 

Given the fact that the number involved in the statistics is small, and the fact that the 
calculations are rather straight forward and do not involve any unknown quantities or quan- 
tities with large uncertainties, the observed and the expected probabilities are certainly 
not inconsistent with each other. 

4.4. Optical depth to microlensing by objects other than the LMC stars 

The probability of microlensing due to stars in the disk and halo of our own galaxy has been 
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discussed by Gould et al. (1994) and the possible contribution of LMC halo was discussed 
by Gould (1992). The probability of the events being caused by the expected number 
of white dwarfs in the halo also appears to be too small. Apart from other theoretical 
reasons, the mass of the objects derived from the observed events do not agree well with 
this being caused by the the white dwarfs in the halo since they have a mass distribution 
which is highly peaked around 0.6 M0 (Weideman 1990). 



5. Can we observationally distinguish between a MACHO and a star in the LMC? 

This probability analysis may be enough for most astronomers to abandon the idea of 
MACHOs since the microlensing can be effectively caused by the known stellar population 
in the LMC itself and one does not need to resort to unknown objects such as MACHOs 
for an explanation. However, let us go one step further and investigate whether there is a 
distinguishing feature which can tell us whether the microlensing is caused by a MACHO 
or a star in the LMC. We will now see that under certain conditions, the fit discrepancy 
near the peak of the light curve can indeed be a distinguishing feature of the microlensing 
being caused by an object within the LMC. Let us see this in detail. 

The amplification sharply rises when the lensing object comes close to being perfectly 
aligned with the source. The physical reason for this is that when the lensing object 
and the source are perfectly aligned, the image becomes a ring (also called the Einstein 
ring) instead of two separate images and thus the amplification rises sharply. [For details, 
see Paczynski, 1986]. In the case of lensing by MACHOs, Re is of the order of 1.2 x 
10^^ \J cm, assuming = 10 kpc and Dg = 55 kpc which, at the source plane, is 

'^10''"^^/-^ cm. The diameter of the source is thus 2 to 3 orders of magnitude smaller 
even if the lensing object is ^O.IM©. Thus the probability of the source and the lens being 
perfectly aligned is extremely small, and can happen only for events with extremely small 
u (;$ 0.01), which is unlikely and certainly not the case in the observed events. However, 
the value oi Re projected to the source plane can be much smaller if the lensing object is 
close to the source. For example, if we consider a typical distance of about 100 pc between 



the source and the lensd, ~ 10^^^ cm, which, for a 0.1 Mq lens, is only about 
about 3 times smaller than the typical radius of a red giant. Thus, in order that at least 
a part of the source is perfectly aligned with the lensing object, the impact parameter has 
to be ^ 0.3, which may explain the fit discrepancy near the peak of the light curve. 

Note also that for the maximum amplification to be appreciable, the source size must not 
be too large for a different reason. If it is too large, even if both objects are perfectly 
aligned, most of the source is still not perfectly aligned with the source, this reduces the 
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maximum amplification. This is given by (Paczyhski 1986) 



Ama.= [l + ^]'^ (12) 

5.1 Light curve due to an extended source 

Assuming the source to be extended, let us calculate the amplifications. A general method- 
ology for such a case was developed by Bontz (1979). The geometry of our particular case 
however (shown in Fig. 1) enables us to develop a simpler solution, which is found to be 
more suitable for numerical integrations. A mathematical approximation and the resulting 
analytical solution for such a case has also been given by Schneider et al. (1992) which, un- 
fortunately for our purpose, does not have sufficient accuracy just at the transition point 
where the fit discrepancy begins to occur. Nevertheless, this analytical expression was 
found to be very useful in providing a constant check on the correctness of the numerical 
algorithm used here to calculate the light curve. 

In a general case, the amplification caused by an extended source would be given by (Eq. 
6.81 of Schneider et al. 1992) 

J(Py I{y) fipjy) 

where I(y) is the surface brightness profile of the source, iJ,p{y) is the amplification of a 
point source at point y, and the integration is carried out over the entire surface of the 
source. 

Let us assume the source to be a disk of uniform brightness with radius ro and let u be 
the impact parameter {=1/Re)- yo be the distance between the lens and the center of the 
source projected onto the lens plane. Let us choose a circular coordinate system and let 
(r, 9) be the representative point which is at a distance y from the lens. From Fig. 1 we 
see that 



y ^ \J Vo + r"^ - 2 yo r cosi'j) (14) 



where 



I.e. 



yo — y Xq + l'^ , ^ = TT — 9 — a and a — sin ^(l/yo) (15) 



y=^yo+r'^ + 2yor cos{a + 9) (16) 
The amplification can be written as 

= ^(^2 +4)1/2 (17) 
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Thus from eq.(13) the amphfication can be expressed as 




(18) 



This is an expression suitable for numerical integration and can be easily integrated. The 
light curve can be reproduced as a function of time t, which is related to Xq through the 
relation xq = {RE/to)t where to is the timescale of microlensing (which is the time taken 
by the lens to cross its own Einstein ring), y in turn, is related to xq through eq.l5 and 



Since I myself do not have the original data, I digitised the published light curve, in order 
to facilitate my attempts on fitting the light curve. Hence the error bars are not shown 
here, although the error bars at the points near the peak were also carefully digitised. The 
light curve was numerically calculated for different values of u and tq to reproduce the 
observed curve. The best fit to the light curve is obtained for u = 0.17 and the value of 
To/ Re between 0.17 and 0.175, which is shown in Fig. 2. The solid curve is the best fit 
obtained using the lens to be a MACHO, and the dotted curve is obtained assuming the 
lens to be a star within the LMC. As is clear, the curve due to a MACHO has a significant 
discrepancy at the peak of the light curve both in the blue and in the red which is reduced 
in case of the best-fit light curve with a star in the LMC as the lens. 

For a analysis, let us first consider the red light curve, the observations of which have a 
higher signal to noise ratio. The total gain, defined as 'E{[A{obs) — A{MACHO)]/a}'^ — 
{[A{obs) — A{LMCstar)]/a}'^ where A refers to the amlification, is about 7.3 units if we 
consider the 4 points near the peak. For the blue light curve, which has a lower S/N, 
the total gain for these 4 points is about 4 units. For other points, both curves are 
practically identical. 

It must be emphasized here that the source is not always extended for an LMC lens. The 
source would be point like for smaller values of vq/Re (which would be the case if the mass 
of the lens is slightly higher, or, the distance between the source and the lens is slightly 
higher). As shown in the next section, the solution becomes fairly unique if the light curve 
can be reproduced by an extended source approximation. 



If this extended source approximation is valid, we can use the value of u and tq/Re as 
calculated in the preceeding section and proceed to calculate some physical parameters. 
The source, from its color and luminosity, has been inferred to be a clump giant, for which 
we can assume the radius to be tq = 7 x 10^^ cm. Now, 
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6. Mass of the lenses and further confirmation of the scenario 



ro ro c 
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I.e. 



-0.3 (20) 



The mass the lens can also be expressed as 



16GD(\ -V?) ^ ' 



From Eq. 20 and 21, which now provide a consistency check, we get 

M D „ . . / K 



= 0.44 " ~ 0.3 (22) 



Mq lOOpc \4:0kms- 

Thus, irrespective of the value of M and D for this particular case, we get the value of Vg 
to be ~ 27 km s~^. The fact that this is similar to the value of dispersion obtained for 
the bar from radio observations (Rohlfs et al. 1984) further confirms the validity of this 
scenario. The mass of the lens can thus be expressed as 

M^O.li^Me (22) 

Although the signal to noise ratio in the EROS events might not allow us to do such a 
detailed analysis of the light curve, the observed amplifications would imply that the source 
can be approximated as a point source. For the point source approximation to be valid, we 
need tq/Re < 1 (Fig. 1) which provides a lower limit to the mass of the lens. An upper 
limit to the mass of the lens comes from equation 21 (using the velocity dispersion). The 
mass range thus derived is shown in Table 1. 

7. Conclusions 

No dark matter is required to explain the observed microlensing events of the LMC stars. 
The probability of microlensing being caused by a star in the bar of the LMC is found to 
be consistent with the observations. Detailed analysis of the light curve of the event with 
the best signal to noise shows that the light curve is better reproduced by a star in the 
LMC than by a MACHO, and indicates that a low mass star in the LMC must be the 
cause of the microlensing. 

The LMC induced events should be strongly clustered towards the central region of the 
LMC. (They should be proportional to the number of monitored stars multiplied by the 
integrated stellar mass density along the line of sight.) In case of the the galactic events, 
for a given number of monitored stars, the events should be uniformly distributed over 
the whole of LMC. (They should be simply proportional to the number of monitored 
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stars.) Furthermore, under certain conditions, the observed hght curve in case of the LMC 
event can be different from the hght curve due to a galactic lens. Thus the two scenarios 
can be distinguished statistically in most cases, and individually in a few cases. Such a 
distinguishing feature seems to have already been seen in case of the light curve of the 
MACHO event. 

Harvey Butcher, Paul Hodge and Bengt Westerlund promptly replied to my queries which 
were helpful in clarifying some doubts regarding the structure of the LMC. The author 
is grateful to Bohdan Paczyhski for his very encouraging and useful remarks. The help 
of Enrique Perez in the digitisation of the light curve was valuable. Some details of the 
observational program provided by Ken Freeman are gratefully acknowledged. 
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Figure captions 

Fig. 1 Schematic of the geometry of the microlensing used for calculation of the light curve. 
The lens is assumed to be a point lens. Note that the Einstein ring radius associated with 
an LMC lens is much smaller than either a disk or a halo lens. Consequently, the minimum 
impact parameter u {= I/Re) needed for a part of the source to be perfectly aligned with 
the lens can be higher. This implies that the source can be resolved at much smaller 
magnifications if the lens resides in the LMC. See text for the details of the symbols used. 

Fig. 2. The observed microlensing events both in blue and red. The best fit due to a 
MACHO is shown by a solid curve. The best fit due to a star within the LMC as the lens 
is shown as the dotted curve. (See text for details.) 



TABLE 1 

Details of the microlensed stars and the derived masses of the lenses 



Event 


Magnitude 


AmplificatiorAmplificatioiEvent 


Impact 


Mass 






(in blue) (in red) 


duration (d) parameter 


(Mo) 


MACHO 


19.6 




33.9 


0.17 


0.1 


EROS 1 


19.3 


liO.lmag liO.lmag 


44 


0.425 


0.002-0.32 


EROS 2 


19.3 


l.l±0.2mag 1.3±0.2mag 53 


~0.34 


0.0002-0.44 



The masses are derived under the assumption that the lensing is caused by stars within 
the LMC. Note that the lower limit to the mass of the lens comes from the requirement 
of a point source as explained in section 6. This is only a lower limit and consequently, 
the actual mass can be much higher than this limit. On the contrary, the upper limit 
comes from the fact that the tangential velocity of the lens cannot be much higher than 
the velocity dispersion (Eq 21). This is an exact expression and considering the usual 
distribution of the velocity dispersion, the mass is more likely to be closer to the upper 
limit. For the first event of Aubourg et al. we have used the radius appropriate for a 
subgiant and for the second event the radius appropriate for a normal star as indicated 
by their colours and luminosities. For a distance D between the source and the lens, the 
mass has to be multiplied by lOOpc/D. For the MACHO event, the amplifications and the 
impact parameter are from the extended source approximation. 
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